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ABSTRACT 



An experimental study was performed to determine the feasibility of 
a laboratory experiment to test an existing theoretical model describing 
sound propagation into a fast bottom underlying a wedge shaped medium. 
Sand under fresh water was found to satisfy the constraints of the theor- 
etical model and to simulate the continental shelf. In a laboratory ex- 
periment, accuracy of density, sound speeds, and attenuation was shown 
to be sufficient to allow quantitative comparison to the predictions of 
the beam angle. 
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I. INTRODUCTION 



The behavior of radiation in a wedge-shaped medium overlying a fast- 
bottom medium has been investigated both theoretically and experimentally 
with both optical and acoustic sources. For example, in 1971 Tien and 
Martin [Ref. 1] examined the behavior of a laser beam coupled into a thin, 
tapered, dielectric film deposited on a substrate with a higher refrac- 
tive index. As the light propagated toward the apex of the wedge shaped 
film, perfect reflection was observed until the changing angle of inci- 
dence decreased below the critical angle, then the light was converted 
into radiation in the substrate. 

An acoustic analog to the above optical study was performed by 
Kuznetzov [Ref.’ 2] in 1973. Kuznetzov developed a normal mode theory 
for sound propagation in both a v/edge and underlying half-space substrate, 
where the substrate sound speed was faster than that of the wedge. 
Kuznetzov concluded that: 

1. Sound traveling toward the wedge apex would be total 1 reflected 
until the angle of incidence decreased to the limiting angle of total 
reflection. 

2. Any sound incident at less than the limiting angle of total re- 
flection would be completely refracted into the underlying half-space. 

This total refraction would occur along the wedge/half-space boundary 
from the wedge apex to the point where the limiting angle was first 
achieved. 

3. Acoustic energy in the half-space would be columnated into a 
well-defined beam. The beam's maximum pressure would occur at an angle 
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of depression (measured from the plane of the wedge/half-space interface) 
that would lie between 8 and 28» where 6 is the wedge angle. 

Kuznetzov performed a series of experiments that supported his theory. 

The formation of a well-defined beam within the substrate (for the 
optical case) was subsequently demonstrated by Tien, Smol insky and 
Martin [Ref. 3]. Reference 3 also presented two theoretical results: 

(1) ray-optics predicted refraction into the substrate beyond the cutoff 
distance (the point at which the critical angle was reached) and (2) 
theory predicted refraction into the substrate before the cutoff distance. 

A similar well-defined beam was predicted and observed by Sigelman, 
et al. [Ref. 4j in a water-aluminum system. Maximum pressure occurred 
at 11 degrees below the water-aluminum interface when the wedge angle was 
1.3 degrees. This was well outside the range predicted by Kuznetzov. 

The observed beam was broader than predicted. 

In 1980, Bradshaw [Ref. 5j extended a computer model first developed 
by Kav/amura and loannou [Ref. 6j to describe the behavior of sound in a 
fast bottom underlying a wedge shaped medium. The model predicted the 
formation of well-defined beams in the bottom, and predicted the effect 
of attenuation in the bottom on the beam. A comparison of the predic- 
tions of this model with the experimental data of Netzorg [Ref. 7] gave 
qualitative agreement, but since the experiment did not match all model 
constraints, quantitative comparison was impossible. It is the purpose 
of this laboratory analysis to determine if a laboratory experiment 
that fulfills the constraints of the model can be designed. If at all 
possible, it is further desired to model the continental shelf in the 
design of the experiment. 
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II. CONSTRAINTS 



The general design of the proposed experiment is illustrated in 
Figure 1, Constraints will be imposed by both the theoretical model 
and the measurements. 

A. MODEL CONSTRAINTS 

The Bradshaw model [Ref. 5j requires that the sound speed in the 
bottom exceed that in the overlying wedge. Physical parameters required 
to be known are the density and sound speed of the two media, the slope 
of the wedge, and the attenuation of the substrate medium. A further 
constraint, dependent on the substrate sound speed, will be discussed in 
detail later. This model assumes planar waves incident on the wedge- 
bottom interface. 

B. OTHER CONSTRAINTS 

Since it was desired to model a real-world environment, specifically 
the continental shelf, a laboratory set-up consisting of fresh water 
over sand was an obvious first choice. The sound speed of fresh water 
was expected to be about 1481 m/s [Ref. 8] and the speed in sand about 
1700 m/s [Ref. 9]. Thus, the requirement of a faster bottom medium 
would be readily fulfilled. 

A wedge of water would require a high frequency (i.e., greater than 
100 kHz) to permit planar incident waves in a laboratory-sized experi- 
ment. If one hundred wavelengths are desired to ensure plane v/aves, at 
100 kHz the source would need to be 1.5 m away from the wedge apex. 
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The use of a water/sand system would permit the hydrophone to be 
moved about in both media; the use of plaster to simulate rock, for 
example, would not. While it would also be necessary to maintain a flat 
slope of at least five degrees between the two media, a water-sand sys- 
tem was expected to be adequately stable. The smoothness of the slope 
is required by the model; the five degree slope encompasses the maximum 
for the continental shelf. 
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Fio-ure 1. Model Geometry 

o 
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III. EXPERIMENTAL DESIGN 



A. MATERIAL SELECTED 

Fresh (tap) water and #30 fine sand were the media used in the experiment. 
The grain size of #30 fine sand varies from 0.70 mm to 0.15 mm. It was 
readily demonstrated that the sand could maintain a slope as large as 25 
degrees. To remove all air bubbles, the water/sand system was heated 
with immersion heaters. Bleach was added at a ratio of one-half gallon 
bleach to 70 gallons of water plus sand in order to control biologic 
growth. 

B. VELOCITY MEASUREMENT TECHNIQUES 

An excellent summary of methods for measuring sound speeds in rocks 
and minerals is given by Anderson and Liebermann [Ref. lOj . Specialized 
methods for use on small specimens, iron compounds, etc. or for high 
pressure and temperature systems were eliminated as inappropriate for 
the experiment being considered. Ultrasonic interferometry was also dis- 
carded despite its exceptional accuracy (0.01 to 0.04 percent) as it is 
best applied to nonporous samples and utilizes frequencies (10-60 MHz) 
that would cause unacceptable attenuation in the verification experi- 
ments. Resonance methods with appropriate frequencies and accuracies 
were undesirable as they are applicable to crystalline specimens only. 
Ultimately, pulse transmissions were selected as they best fit the 
necessary criteria. This method is applicable to either fluid or por- 
ous materials, acceptable over a frequency range of 50 kHz to 10 MHz, 
and accurate to within one to three percent. 
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1. size Criteria 



Steel-bound glass tanks measuring 70 cm square by 60 cm high 
were readily available. Length vs depth criteria, from Reference 10, 
require that the layer of sand on the tank bottom should not be less 
than 14 cm to avoid boundary reflection problems. To satisfy this con- 
straint, a sand layer of 20 cm was chosen. 

To avoid pulse dispersion. Reference 10 also suggested that the 
sand depth not be less than five wavelengths. Since the sound speed of 
the sand was expected to be 1700 m/s, and the frequency would be of the 
order of 100 kHz, five wavelengths would be 8.5 cm. A sand depth of 
20 cm would also satisfy this restraint. 

To avoid particulate scattering, the wavelength must be at 
least three times the grain size [Ref. 10]. Since the largest grain 
size in number 30 sand is 0.070 cm, a wavelength, equal to 1.7 cm, 
readily qualifies. 

Since sound speed is not a function of frequency [Ref. 9], re- 
sults from this laboratory analysis were expected to be applicable to 
frequencies used in the experiments. 

2. Electronic Equipment 

A schematic of the equipment configuration is shown in Figure 2. 
All components were off-the-shelf and readily available. 

Output from a General Radio model 1310 oscillator with a fre- 
quency range of 2 Hz to 2 MHz was fed simultaneously into a frequency 
counter and a tone burst generator. The counter, a Hewlett-Packard 
5233L, would read +_10 Hz at 100 kHz. The General Radio Type 396-A tone 
burst generator was used to generate either 8 or 16 cycle pulses. The 
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output passed through a Hewlett-Packard HP 467-A power amplifier before being 
fed to a transducer. The pulsed signal was also fed to the trigger in- 
put of a Tektronix type Rm-503 oscilloscope. 

Three types of transducers were used as sources. The first was 
a homemade, mylar, broad-band transducer with a 3 cm by 9 cm active 
face. This source was quite directional at 100 kHz but proved to be 
difficult to move through the sand (for sound speed and attenuation 
measurements) and was not sufficiently powerful to generate usable sig- 
nals in the sand (as would be required in model verification runs). The 
mylar transducer required by 150 V DC supply and polarizing network. 

A second, smaller, homemade transducer was tried as a source. 

This transducer, a matrix of ceramic elements, has a 1 x 2 cm active 
face and was highly directional at 100 kHz. This transducer offered 
little improvement over the first. 

Celesco Industries type LCIO hydrophones had been consistently 
used as receivers. The LCIO is a small (0.97 cm diameter by 2.87 cm 
long) cylinder, with a receiving range of 0.1 to 120,000 Hertz. The 
transducer is designed to be omnidirectional in a plane perpendicular 
to the axis of the cylinder with a tolerance of +^1 dB at 100 kHz, and 
omnidirectional in a plane containing the axis of the cylinder with a 
tolerance of +2 dB at 25 Hz. The omnidirectional characteristic was 
not desirable for the source, as this would increase the number of pos- 
sible reflected signals from the walls, water surface, etc., but the 
small size and higher signal strength were desirable. Geometry was used 
to isolate desired signals and will be described when appropriate. 
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Because the direct and reflected pulses arrived from different 
directions, use of LClO's as the source and receiver for attenuation 
measurements initially produced non-reproducible results. Directivity 
measurements in the plane perpendicular to the axis of the cylinder were 
made and the results are shown in Figs. 3 and 4. Subsequent measurements 
used LCIO number 2338 as a source and LCIO number 2319 as a receiver. 

The receiver was oriented to use the area about 50 degrees as the receiv- 
ing face. The transducers were always arranged geometrically to use the 
radial plane. 

The received signals v/ere amplified 20 dB by a Hewlett-Packard 
HP-465A amplifier, then passed through a Spencer-Kennedy Laboratories, 
Inc. model 302 variable electronic filter (set at 60 kHz high pass) to 
eliminate low frequency mechanical noise present in the laboratory be- 
fore being passed to the oscilloscope. The data of Fig. 5 demonstrates 
that spherical spreading (1/r) was observed for all LCIO source-receiver 
separations greater than 8.4 cm. All measurements were made under far- 
field conditions. 
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Transducer 
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Figure 2. Electronic Equipment Schematic 
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Figure 3. LCIO 2338 Directivity Pattern 
(units of volts) 
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Figure 4. LC10 2319 Directivity Pattern 
(units of volts) 
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IV. MEASURED DATA 



For all of the following measurements, the sand layer in the tank 
was kept flat and horizontal, as only physical parameters were being de 
termined and not wedge characteristics. 

A. DENSITY 

1 . Water 

According to Lange's "Handbook of Chemistry" [Ref. Ill, the 

density of distilled water ranges from 0.99913 g/cm^ at 15°C to 0.99707 
3 o 

g/cm at ZS°C. The expected density of room temperature water, to 3 
significant figures, was therefore 1.00 g/cm . Measured volumes of the 
bleached water were periodically weighed, and the observed density was 

3 

1.00 g/cm , with no noted variation. 

2. Sand 

The density of water saturated sand was measured by partially 
filling a v/eighed 100 ml graduated cylinder with a water over sand mix- 

3 

ture, observing the volume of each, and using a density of 1.00 g/cm 
for the water to calculate the weight of the water volume. The density 
of the water saturated sand, for eight separate measurements, was 
1.98;^0.03 g/cm^. 

B. SOUND SPEED 
1 . Water 

Three separate sets of sound speed in water measurements were 
made. The first involved the use of the mylar transducer as the source 
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One LCIO was placed 20.1 cm away, and the second LCIO was moved along 
the straight line defined by the source and first receiver. The time 
of flight between the two receivers was measured. A simultaneous check 
for 1/r dependence confirmed these measurements were all in the far- 
field. The data are summarized in Table 1 of Appendix A. The average 
sound speed was 1454 m/s with an (N-1) variance of 33 m/s. 

The second set of data utilized a reflected signal to increase 
the total path length. The mylar source was fixed, and the receiver, 
placed between the source and the wall, was moved along a perpendicular 
line between the source and the wall. The time of flight was the time 
difference between the arrival of the direct and reflected signals; the 
distance was twice the distance from the receiver to the wall. These re- 
sults are summarized in Table 2 of Appendix A. The average speed was 
1459 m/s, with a variance of 14.7 m/s. 

In the third set of measurements, the distance over which the 
time of flight was measured was fixed. The LCIO source was suspended 
11.0 cm from the tank wall, the second LCIO, placed further from the 
wall than the source, was moved along the perpendicular formed by the 
wall and the source. As the receiver was moved, the distance differ- 
ence between the direct path signal and the signal reflected from the 
wall remained a constant 22.0 cm. The time difference between the ar- 
rival of the two signals was recorded. Table 3 of Appendix A summarizes 
the results. The average speed was 1460 m/s, the variance was 19 m/s. 

The data compared well with theory. From Kinsler and Frey 
[Ref. 8], the equation of sound for distilled water is given by 

c = 1403 + 5t - 0.06t^ + 0.0003t^ (1) 

w 
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where t 1s the temperature 1n degrees Celcius. A comparison of measured 
and theoretical values is presented in Table VI-1. Note that direct 
path results showed greater variation, but all results were within two 
percent of expected values. 

2. Sand 

Five separate sets of sound speed measurements were made in the 
sand; four in conjunction with attenuation data, and one to check for a 
sound speed gradient. All three source transducers were used. For the 
attenuation runs, the source was buried at a fixed depth, probing was 
with an LCIO receiver inside a glass tube for the depth of maximum signal 
amplitude. The glass tubes were necessary to prevent abrasive damage to 
the receiver. The data are summarized in Tables 1 through 4 of Appendix 
B. 

The fifth set of data utilized the two LClO's fixed 30 cm apart. 
The transducers were initially placed in water and sand was added to 
progressively bury them. The entire process was carried out underwater, 
so degassed, saturated sand was used. The data are summarized in Table 
5 of Appendix B. Table IV-2 lists the resultant sound speeds. 

a. Theoretical Sound Speed in Sand 

Since a search of the literature revealed many theoretical 
and experimental values for the speed of sound in sand under brine, but 
not under fresh water, and since the values given were greater than 
those observed in this experiment, the following analysis was performed 
to derive a theoretical value for sound speed in sand under fresh water. 

From Urick [Ref. 12], the density of a 1 iquid -sediment 
mixture is given by 
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Table IV-1. Measured vs Theoretical Sound Speed in Water 
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( 2 ) 



Omix ' (' - 

where = the density of the mixture, 

Y = the porosity, 

= density of water. 



and p^ - the density of individual sand grains. 

From equation (2), using the measured value p . =1.98 g/cm^ and 

mi X ^ 

3 

p^ = 2.69 g/cm from Reference 9, gives y = 0.42. Also from Reference 
12, the speed of sound in the saturated sediment is given by 

"^mix " Ps^^ • (3) 

where and are the compressibility of the water and sand, respec- 
tively. Values for k^ and k^ were calculated from the relationships 



and 

Setting c, , 
^ w 

k. , yields 
w 



*"w ^*^w^w^ 

S = (Ps^) 



• 1/2 



• 1/2 



= 1481 m/s [Ref. 8], and = 1.0 g/cm , and solving for 



(4) 

(5) 



k^^ = 4.56 X 10"^^ m sec^ g"^ (6) 

The sand in the experiment was composed of quartz particles. From 
Anderson and Lieberman [Ref. lOj, page 360, c^ was set to 7000 m/s. 
Solving equation (5) yields 

kg = 7.59 X 10"^^ m sec^ g"^ (7) 

Substituting the above values in equation (3) and solving for c^^^ 

yields an expected speed for sand under fresh water of 1605 m/s. 

A comparison of the expected value with the sound speeds 
summarized in Table IV-2 shows close correlation with the value obtained 
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TABLE IV-2. Measured Sound Speeds in Sand 



Measurement 

Method 


Sound Speed 
(m/s) 


(N-1 ) a 
(m/s) 


Variable 


distance 


1579 


58.6 


Variable 


distance 


1597 


44.7 


Variable 


distance 


1434 


71 .9 


Variable 


di stance 


1567 


32.6 


Fixed source/receiver 


1613 


0.0 
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using a fixed source-receiver distance (1613 m/s) but not close agreement 
with the majority of values obtained from measurements utilizing glass 
tubes. The use of glass tubes disturbs the medium, and the distance from 
source to receiver is not accurately measurable. Over the short distances 
used, a 0.5 cm error in the distance would lead to significant errors. 

Thus, the value obtained from the fixed source and receiver method is more 
likely accurate, and, not unexpectedly, shows closer agreement with the 
theoretical values. 

b. Sound Speed Gradient 

As illustrated in Fig. 6, Hamilton [Ref. 13] predicts a pro- 
nounced sound speed gradient in a fine sand sediment. A change in slope 
of 27 s“^ is specified at zero depth. Although the gradients were estab- 
lished for sand under brine, they are attributed to a decrease in sediment 
porosity with pressure, to temperature effects, and to pressure effects on 
the pore water, and would therefore be expected to be present in fresh 
water systems as well. Sound speed measurements tabulated in Table 5 of 
Appendix B were taken specifically to check for the presence of a gradient. 
None was found within the accuracy limitations of the method. For a sound 
speed of 1613 m/s, a 2?o error (see Sec. III-A) would be 32 m/s and would 
mask a 27 m/s change. Even if a 27 m/s gradient were present, the re- 
sulting 62 m radius of curvature would not be sufficient to be detectable 
over the experimental ranges. The lack of a detectable gradient, although 
not mirroring actual continental shelf conditions, is consistent with the 
isovelocity substrate assumed in Bradshaw's model [Ref. 5]. 
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Velocity, m/s 




Figure 5 . Sound Speed Gradient in Brine 
Saturated Sand (adapted from 
reference 13.) 
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C. ATTENUATION 



1 . Model Limitations 

It was mentioned in Section II-A that a further constraint, to be 
defined once the substrate sound speed had been determined, must be ful- 
filled for verification of Bradshaw's model. This constraint is that the 
signal must propagate beyond a specific range determined as follows. The 
model uses a Green’s function analysis to obtain the pressure and phase 
distribution of the sound in the bottom. To solve the appropriate Green's 
function, it was necessary to assume that 

!< 2 lr - r'\ » 2tt (8) 

where l <2 = u/c^ is the wave number for the sound in the sand. Using the 
determined sound speed of 1613 m/s and assuming a frequency of 100 kHz, 
k 2 equals 389 m"\ Thus, it is required that |r - r"] » 0.016 m 

The quantity |r - r'\ is defined in Fig. 7. The quantity r' is 
the distance from the wedge apex to a point on the wedge-substrate 
boundary (integration in the model is carried out over r'). The quantity 
f is the distance from a point along the wedge-substrate boundary, called 
X, to the field point. The distance X is the distance at which the cri- 
tical angle is reached and where lowest mode transmission into the sub- 
strate begins. The distance is defined by 

^ 4 f sin9^ sing’ 

9^ = arccos (c^c^) 

where f is the frequency, 

6 is the wedge angle 

Cw» c^ are the sound speeds in the wedge fluid and substrate, 
respectively. 



(9) 

( 10 ) 
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Figure 7. Geometry ror the Wedge Model. 



and 9^ is the critical angle. 

For f = 100 kHz, a wedge angle of 3% c , = 1469 m/s and c^ = 1613 m/s, X 

w s 

would equal 0.17 m. 

Typically, within the model, the integration limit is from nX to 
NX, where n = 0,1,2,..., N = 1,2,3,..., N>n. In the simplest case, the 
integration would be carried out from 0 to X. This places a constraint 
on r, and thus on permitted attenuation. 

Consider, as a minimum, a factor of 10 to satisfy the "much 
greater than" inequality; this would require |r - r'\ to be greater than 
0.16 m. This implies mathematically that no measurements may be made at 
point r within a circle of radius 16 cm about the point r", and since r' 
can vary from 0.0 to 0.17 m, an area of "forbidden territory" is traced 
about the wedge apex and boundary surface. Thus, to test the model, 
since measurements must be made outside this area, the attenuation must 
not be so great that insufficient signal remains beyond the forbidden zone. 
(If "much greater than" were taken to be 100 times, measurements would 
have to be made over 1.6 m from the boundary). 

2. Results 

The methodology for collecting attenuation data was described in 
Section IV-B-2. The data were processed graphically. Taking the 
natural logarithm of the well-known equation 

V- -ar 

V = -°e (11) 

r 

where V is the measured voltage, 
is the source voltage, 
r is the distance from the source, 
and a is the attenuation. 
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the linear relation 



In(Vr) = In(V^) - ar (12) 

is obtained. Thus a, in nepers per meter, will be the slope of a graph 
of In(Vr) vs r. Graphs of the four data sets of Appendix B are shown in 
Figures 8 through 11. A summary of the attenuation values is presented 
in Table IV-3. 

Figure 12 shows a plot of the attenuation vs frequency. The line 
indicates the slope of any line having a dependence of attenuation on the 
first power of frequency. All data fall within the limits of data col- 
lected by Hamilton [Ref. 14j. 

3. Reflection Coefficients and Signal Level 

Normal incidence pressure reflection and transmission coefficients 
were measured as an accuracy check on the measured p and c values, and 
also to provide an estimate of the sound pressure level that would be 
transmitted from the wedge to the underlying medium. A summary of the 
measured data is presented in Appendix C. 

The reflection coefficient R is calculated from 



R 



= iPw^w- PsS 
iPw^w-" PsSi 



and the transmission coefficient T from 



(13) 



T = 1 - R (14) 

With = 1.00 g/cm^, = 1.98 g/cm^, c^ = 1469 m/s, and c^ = 1613 m/s, 
R is calculated to be 0.37 and T to be 0.63. The measured values of re- 
flection were 0.31 and 0.29. The measured transmission coefficient was 
0.65. These values indicate that the measured density and sound speed 
values are within the expected range. 
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Figure 8. Attenuation in Sand 
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Figure 9. Attenuation in Sand 
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Figure 10. Attenuation in Sand 
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Figure 11. Attenuation in Sand 
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Figure 12. Attenuation vs Frequency 
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From Table 1 of Appendix C, a received voltage of 6 volts is 
typical. Assuming a desired reading on the order of 1 volt for test 



runs, a drop in level of no more than 15.6 dB is indicated, 
the measured attenuation of 56.6 dB per meter would indicate 
ments could be made at least 27 cm from the wedge boundary, 
fies the limit imposed on r in Section IV-C-1. 



At 100 kHz, 
measure- 
This satis- 
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TABLE IV-3. Attenuation Results 



f(kHz) 


a(Ne/m) 


a(dB/m) 


Correlation 

Coefficient 


75.10 


7.26 


63.2 


-0.90 


100.00 


6.50 


66.6 


-0.99 


100.04 


7.39 


64.3 


-0.76 


130.00 


9.79 


85.2 


-0.92 
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V. ERROR ANALYSIS 



It is important to know if the physical properties were determined 
with sufficient accuracy to verify Bradshaw's model. Bradshaw developed 
an analytical equation for the depression angle of the beam formed in 
the substrate. As this analytical equation was based on model results, 
and since the equation is dependent on all measured physical properties, 
the equation should provide a realistic basis for error analysis. The 
analytical equation is 



6p = 17.22 



gO.329 ,^0.772 



(15) 



where Gq is the beam depression angle in degrees, 

S is the wedge angle in degrees, 

0^ is the critical angle in radians, 
and Pw/^s ratio of the density of water to the density of 

saturated sand. 

The critical angle is defined by 

0^ = arccos (c, ,/c.) (16) 

c w' b 

where c and c are the sound speeds in water and sand, respectively, 
w s 

Using standard error analysis, and noting p^/p^ = 1/p^, 



( 17 ) 



= 5.66 0p (I/P 5 ) 

+ 13.29b‘^-^^^ (I/P 5 ) 



0.25430.671^3 



0.254„ 0.228.. 
®c "‘^c 
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Similarly 



+ 4.378°-32S e“-”2o 

c s s 



d0 = dc + dc 
*- TF~ w -7-r- s 
“'w 



(18) 



dc, - 

H7 ' , 

s s 



1 dc 1 

^ 1 



^ - (VS)' 



- 1/2 



where c^ - 1469il5 m/s from Table 2, An error of 3 percent, in keeping 

with the worst case predicted accuracy from Reference 10, was assigned to 

sound speed in saturated sand. Therefore, c = 1613^48 m/s and dG = 0.04 

s c 

radians . 

Assuming the wedge angle can be measured to within 0.5 degrees, the 
following values were used to evaluate equation 17. 

= 1.98 g/cm^ 

S = 3° 



dp^ = 0.03 g/cm" 



0^ = 0.4257 rad 



dS = 0.5 



d0^ = 0.04 rad 



Therefore 



d0Q = 3.64 + 0.75 + 0.16 
= 4.55*’ 

0Q for the above parameters is 15.21 degrees. The deviation of 0.75° 
attributable to sound speed measurements is roughly 5 percent, but is 

o 

dwarfed by the 3.64 (24 percent) error attributable to the wedge angle 

measurement. A method of accurately measuring and controlling the 
wedge angle will be a necessity for successful model test runs. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A laboratory sized experiment to test the predictions of Bradshaw's 
model is possible, using off-the-shelf equipment in a water-over-sand 
system. Measurements can be made of all necessary physical parameters, 
and attenuation is not so severe as to preclude measurements at suffi- 
cient distances to satisfy model constraints. 

Densities and sound speeds were obtained that agreed with literature 

or theoretical values. The densities of water and saturated sand were 
3 3 

1.00 g/cm and 1.98 g/cm , respectively. The sound speed in water was 
1469 m/s; the sound speed in saturated sand was 1613 m/s. 

The use of a rigid array of receivers is highly recommended for all 
measurements within the sand in order to control source-recei ver spacing. 

Model results are strongly dependent on the wedge angle. The sand 
will maintain a flat, sloped surface, but great care in the formation and 
measurement of the wedge angle will be required. 

It was found necessary to maintain a high bleach concentration in 
order to control gas producing biologic growth. Every change of water 
should be accompanied by a renewal of the bleach level. 
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APPENDIX A 



SOUND SPEED IN WATER - DATA 

TABLE 1. Two Receivers Direct Path Data 

Source Transducer = Homemade Mylar 



Separation 

(cm) 


Time 

(10“^ sec) 


Sound Speed 
(meters/second) 


6.8 


4.50 


1511 


11.7 


7.95 


1472 


16.7 


11.70 


1427 


19.0 


13.30 


1428 


21.4 


14.90 


1436 


26.3 


18.10 


1453 



c = 1454 meters/second 
a = 32.6 (2.2%) 
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TABLE 2. Reflected Path, Varied Distance - Data 
Source Transducer = Homemade Mylar 



Separation 

(cm) 


Time 

(10‘^ sec) 


Sound Speed 
(meters/second) 


28.8 


19.55 


1473 


38.8 


26.32 


1474 


48.8 


33.25 


1468 


58.8 


40.15 


1464 


68.8 


47.85 


1438 


78.8 


53.70 


1467 


88.8 


59.60 


1490 


98.8 


67.00 


1475 



c = 1469 meters/second 
a = 14.7 (1.0%) 
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TABLE 3. Reflected Path, Fixed Distance - Data 



Source Transducer = LCIO #2338 



Separation 

(cm) 


Time 

(10'^ sec) 


Sound Speed 
(meters/second) 


22.0 


14.92 


1474 


22.0 


15.32 


1436 


22.0 


14.88 


1478 


22.0 


14.92 


1474 


22.0 


14.92 


1474 


22.0 


14.92 


1474 


22.0 


15.32 


1436 


22.0 


15.32 


1436 


22.0 


15.08 


1459 



c = 1460 meters/second 
a = 18.8 (1.3%) 
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APPENDIX B 

SOUND SPEED AND ATTENUATION IN SAND - DATA 



TABLE 1. Sound Speed and Attenuation in Sand 



Source Transducer = Homemade Mylar 

Frequency = 100.04 kHz 

Source placed 10 cm below sand surface. 



Distance 




Time 


Amplitude 


Sound Speed 


(cm) 




(10~^ sec) 


(mV) 






(meters/second) 


15.15 




9.04 


11.0 






1676 


16.15 




10.10 


19.0 


(18.0)* 


1599 


19.15 




12.25 


7.0 


( 


7.2) 


1563 


19.60 




12.25 


9.0 






1600 


22.05 




13.80 


7.0 


( 


7.2) 


1598 


23.75 




15.40 


3.0 


( 


2.8) 


1542 


24.10 




17.55 


7.5 






1373** 


25.55 




15.95 


4.0 


( 


4.0) 


1602 


27.65 




17.00 


3.5 


( 


3.2) 


1626 


28.15 




19.15 


3.0 


( 


2.4) 


1470 


29.10 




18.05 


3.0 






1612 


29.70 




18.10 


6.0 






1641 


31.90 




20.20 


1.0 






1579 


31.90 




21.25 


1.5 






1501 


34.20 




22.85 


1.5 






1497 


38.00 




no signal 


— 








*changed 


scale 




c = 


1579 meters/second 


**omitted 

(outside 


point 

2a) 




a = 


58.6 (3.7%) 
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TABLE 2. Sound Speed and Attenuation in Sand 



Source Transducer 


= Homemade Mylar 






Frequency = 100.00 


kHz 






Source placed 14.5 


cm below surface 


of sand. 




Distance 


Time 


Amplitude 


Sound Speed 


(cm) 


(10“^ sec) 


(mV) 


(meters/second) 


10.30 


6.38 


27.0 


1614 


18.39 


11.60 


10.0 


1585 


19.15 


12.12 


42.0* 


1580 


23.71 


14.26 


4.5 


1663 


29.22 


18.94 


2.8 


1543 


*point omitted 
for attenuation 




c = 1597 meters/second 
a = 44.7 (2.8%) 
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TABLE 3. Sound Speed and Attenuation in Sand 



Source Transducer = Homemade Ceramic 
Frequency = 75,1 kHz 

Source placed 13.0 cm below sand surface. 



Distance 

(cm) 


Time 

(10'^ sec) 


Ampl i tude 
(mV) 


Sound Speed 
(meters/second) 


12.2 


9.05 


25.0 


1348 


12.5 


8.00 


20.0 


1562 


23.5 


16.50 


8.0 


1424 


20.6 


14.35 


3.0 


1436 


27.4 


19.70 


2.5 


1391 


32.3 


22.35 


2.5 


1445 



c = 1434 meters/second 
a = 71.9 (5.0%) 
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TABLE 4. Sound Speed and Attenuation in Sand 



Source Transducer = Homemade Mylar 

Frequency = 130.00 kHz 

Source placed 7.5 cm below sand surface. 



Distance 

(cm) 


Time 

(10"^ sec) 


Ampl i tude 
(mV) 


Sound Speed 
(meters/second) 


12.0 


7.7 


10.0 


1558 


13.0 


8.0 


— 


1625 


18.7 


12.0 


3.0 


1565 


23.0 


14.9 


1.8 


1544 


26.1 


16.5 


0.8 


1582 


27.3 


18.1 


0.6 


1508 


31.5 


20.2 


0.6 


1559 


32.2 


20.2 


0.6 


1594 


32.6 


20.8 


0.4 


1571 



c = 1567 meters/second 
a = 32.6 (2.1%) 
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TABLE 5. Sound Speed Gradient Data 



Source Transducer = LCIO #2338 
Frequency = 120.05 kHz 

Source and receiver fixed at 30.0 cm separation. 



Hater 

(cm) 


Sand 

(cm) 


Time 

(10'^ sec) 


Amplitude 

(V) 


Sound Speed 
(meters/second) 


n.o 


— 


20.2 


5.4 


1485 


27.0 


— 


20.2 


10.4 


1485 


30.0 


0.0 


18.6 


0.5 


1613 


29.0 


1.0 


18.6 


1.0 


1613 


27.0 


3.0 


18.6 


2.5 


1613 


25.5 


4.5 


18.6 


1.35 


1613 


24.0 


6.0 


18.6 


1.35 


1613 


23.0 


7.0 


18.6 


1.35 


1613 


21.0 


9.0 


18.6 


1.35 


1613 


20.0 


10.0 


18.6 


1.35 


1613 


17.5 


12.5 


18.6 


1.35 


1613 


15.0 


15.0 


18.6 


1.35 


1613 


12.5 


17.5 


18.6 


1.35 


1613 


10.0 


20.0 


18.6 


1.35 


1613 
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APPENDIX C 



REFLECTION AND TRANSMISSION COEFFICIENTS - DATA 
TABLE I. Surface Reflection Coefficient 



Di rect 


Path 


Reflected 


Path 


Reflection Coeff 


Distance 


Amp. 


Di stance 


Amp. 




(cm) 


(V) 


(cm) 


(V) 




22.5 


5.8 


48.5 


2.9 


1.077 


30.0 


5.4 


37.2 


4.0 


0.9185 


20.0 


6.2 


40.0 


3.3* 


1.010 


20.0 


6.4 


40.0 


3.0** 


0.9880 



R = 0.998 
a = 0.06 

* Directivity (.241/. 254) 

** Directivity (.254/. 241) 

* and ** measured at normal incidence 
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TABLE 2. Bottom Reflection Coefficient, Normal Incidence 



Distance 

Direct = 20.0 cm 
Reflected = 40.0 cm 



Amplitude (V) 



Di rect 


Reflected 


5.4 


0.8 


5.7 


0.8 


6.5 


0.8 


7.2 


1.0 


7.1 


1.0 


7.1 


0.9 


7.0 


1.0 


7.0 


0.9 


7.2 


1.2 


10.0 


1.2 


8.6 


1.1 



Directivity Index 
Direct = 0.254 
Reflected = 0.241 

Refl ection 
Coefficient 

0.31 
0.30 
0.26 
0.29 
0.30 
0.27 
0.30 
0.27 
0.35 
0.25 
0.27 

R = 0.29 
a = 0.028 
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TABLE 3. Bottom Reflection Coefficient, Normal Incidence 



Distance 

Direct = 19.2 cm 
Reflected = 40.4 cm 



Directivity Index 
Direct = 0.254 
Reflected = 0.241 



Amplitude (V) 



Reflection 

Coefficient 



rect 


Reflected 




7.2 


1.1 


0.34 


7.0 


1.2 


0.38 


7.5 


1.25 


0.37 


7.0 


1.2 


0.38 


7.1 


1.0 


0.31 


7.6 


1.0 


0.29 


7.6 


1.0 


0.29 


7.5 


0.8 


0.23 


6.5 


0.8 


0.27 


5.6 


0.8 


0.32 


6.6 


0.9 


0.30 


7.4 


0.8 


0.24 


7.4 


0.95 


0.28 


7.5 


1.0 


0.30 


7.0 


1.0 


0.32 


7.0 


1.1 


0.35 


7.5 


1.2 


0.35 


7.0 


1.1 


0.35 


7.0 


0.9 


0.29 



R = 0.31 
a = 0.044 



53 



TABLE 4, Transmission Coefficient, Normal Incidence 



Direct Path = 20.0 cm in water 

Transmitted = 10.0 cm in water, 10.0 cm in sand 



Frequency 

(kHz) 


Di rect 


Amplitude 

(volts) 

Transmitted 


Transmission 

Coefficient 


125 


8.7 


6.8 


0.78 


130 


9.0 


7.0 


0.78 


125 


8.7 


5.0 


0.57 


120 


7.8 


5.0 


0.68 


115 


7.4 


3.6 


0.49 


110 


6.9 


3.6 


0.52 


120 


7.8 


5.4 


0.69 


125 


8.7 


5.2 


0.60 


130 


9.0 


7.6 


0.84 


125 


8.7 


5.2 


0.60 


120 


7.8 


5.2 


0.67 


115 


7.4 


3.8 


0.51 


115 


7.4 


3.8 


0.51 


120 


7.8 


5.2 


0.67 


125 


8.7 


5.3 


0.61 


130 


9.0 


7.6 


0.84 



T = 0.65 
a = 0.116 
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